QUANTUM CURVE IN g-OSCILLATOR MODEL 



S. SERGEEV 

Abstract. A lattice model of interacting q-oscillators, proposed in [V. Bazhanov, S. 
Sergeev, |arXiv:hep-tli/0509181|, is the quantum mechanical integrable model in 2+1 
dimensional space-time. Its layer-to-laycr transfer-matrix is a polynomial of two spectral 
. parameters, it may be regarded in the terms of quantum groups both as a sum of sl(N) 

. transfer matrices of a chain of length M and as a sum of sl(M) transfer matrices of a chain 

of length N for reducible representations. The aim of this paper is to derive the Bethe 
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Ansatz equations for the q-oscillator model entirely in the framework of 2+1 integrability 
providing the evident rank-size duality. 



Introduction 

The g-oscillator lattice model was formulated recently in [1,2]. It describes a system of 

interacting g-oscillators situated in the vertices of two dimensional lattice, and therefore 
(N ■ 

' it is the quantum mechanical system in 2 -|- 1 dimensional space-time in the same way, 
00 '. 

Tj- ■ as a chain of interacting particles (or spins) is regarded as a model in 1 -|- 1 dimensional 

o ■ 

Q I space-time. Formulation of the g-oscillator model provides a definition of a layer-to-layer 

V/-^ . transfer matrix as a polynomial of two spectral parameters. This transfer matrix may be 

^ ■ interpreted in the terms of quantum inverse scattering method and quantum groups, so 



that both sizes of the two dimensional lattice may be interpreted as either a length of an 
effective chain or as symmetry group's rank. This was called in [1] the "rank-size" duality. 



^ I The appearance of a complete set of fundamental transfer matrices for lAq{sl) series is a 
■ signal that the layer-to-layer transfer matrix of g-oscillator model is closely related to 
Bethe Ansatz in the form of generalized Baxter's "T-Q" equations. The subject of this 
paper is the derivation of such equations in the framework of 2-1-1 dimensional integrability. 

Below in this Introduction we formulate the answer, i.e. we give an explicit form of 
"T-Q" equations in the terms of a given layer-to-layer transfer matrix. To do this, we 
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2 S. SERGEEV 

need to repeat the structure of the layer-to-layer transfer matrix for g-oscillator model in 
more details. 

The g-oscillator model describes a system of interacting g-oscillators ?iv, 

(1) K : x,y, = 1 - , y^x, = 1 - g^hv ^ ^^qh. = ^h.+i^^ ^ y^^h. ^ ^h.-iy^ ^ 

situated in the vertices v of two dimensional square lattice of sizes N x M. Index v 
stands for a coordinate of a vertex. Oscillators from different vertices commute (what is 
called "locality"), the whole algebra of observables is thus 7^*^^^, and the vertex index 
V corresponds to the number of component of the tensor power. In this paper we imply 
mostly the Fock space representation of g-oscillators. 

The two dimensional lattice may be identified with a layer (or a section) of three 
dimensional cubic lattice, further we call it either the layer or the auxiliary lattice. 

Auxiliary matrices La,i3[Hv], acting in ® J^^, were introduced in [1]. The layer 
transfer matrix T{u, v) may be constructed as a trace of an 2d ordered product of auxiliary 
matrices L[Hv]. The transfer matrix is a polynomial of two spectral parameters, 

N M 

(2) T{u,v) = 

n=0 m=0 

its coefficients tn,m € 7^"^^^^ form a complete commutative set. Matrices L[H^] depend on 
some extra C- valued free parameters, for their generic values the model is inhomogeneous. 
The layer transfer matrix Q may be identically rewritten in two ways, 

N M 

(3) T{u,v) ^ Yl = E ^""Tit^'^^) ' 

n=0 m=0 

where 

M 

(4) = E ^'"t--- 

m=0 

is the 2d transfer matrix for h{q{slj^) chain of the length M, corresponding to the fun- 
damental representation tc^^^ in the auxiliary space (here Un stand for the fundamental 
weights of An^i, tt^q and tt^^ are two scalar representations, Tq"*'^-* and T^^^'' may be 
written explicitly). The same layer transfer matrix T{u,v) may be rewritten as the sum 
olUq{slM) transfer matrices 

N 

(5) T(f-)(«) = «"t„,^ 

n=0 

for the length N chain (the last part of (jH))). 
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The result of this paper is the derivation of the dual Bethe Ansatz equations for the 
g-oscillator model. They may be formulated as follows. Let normalized transfer matrices 
be 

(6) r^:i-'\u) = T^y{-i-qru) , tI:^-\v) = T£'")(-(-g)"t;) . 

And let C-numerical parameters of g-oscillator lattice are inhomogeneous enough. 
Then "T-Q" equation for sIm is 

M 

(7) Y.^-vr rl:!^'\u) Qiq'^u) = 0. 

m=0 

The statement is that if tn,m take their eigenvalues, then there exist M special values 
Vi, . . . , vm of V, such that corresponding Qi{u), . . . , Qm{u) in ((7j) are polynomials^. De- 
grees of the polynomials are uniquely defined by certain occupation numbers of oscillators. 
In its turn, equivalent "T-Q" equation for sIn is 

N 

(8) Y.^-ur T^'-\v)Q{q'-v) = 0. 

n=0 

If tn,m take their eigenvalues, then there exist special values ui, ... ,un of u, such that 
corresponding Qi{v), . . . ,Qj^{v) in (jHl) are polynomials. All the other forms of nested 
Bethe Ansatz equations follow from (|7j) or (jHJ. 

Polynomials Q{u) and Q{v) may be denoted in the quantum mechanical way as "wave 
functions" of states {Q\ and \Q): 

(9) Q{u) = {Q\u) and Qiv) = {v\Q) , 
where \u) and (^1 serve the simple Weyl algebra W: u v = v u, 

u\u) = \u)u , v|u) = \q'^u)v , 

(10) and 

{v\u = u{q^v\ , (f |v = v{v\ . 

Let 

N M 

(11) JM = $^$^(-g)-"'"(-ur(-v)-t„,^. 

n=0 m=0 

Then ((7j) and (jHI) are correspondingly 

(12) (Ql J(u,v) \u) = {v\ J(u,v) IQ) = 0. 

"'^parameter v in the "u-shift" equation Q is irrelevant since a re-scaling Q{u) u'^Q{u) changes it. 
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The formulation of the g-oscillator model and definition of T{u,v) are locally 3d in- 
variant, the quantum group interpretation Q is the secondary one. In this paper we will 
derive (fTTj) without any quantum group technique. 

To explain our method, we need to comment a little the classical limit. 

In the classical limit g — ^ 1, the local g-oscillator generators become the classical dy- 
namical variables, the g-oscillator model becomes a model of classical mechanics, quantum 
evolution operators become Baecklund transformations for the dynamical variables. In 
particular, T{u, v) may be understood as a partition function of a completely inhomoge- 
neous free fermion six-vertex model on the square lattice (but it should not be regarded as 
a model of statistical mechanics). In its turn, J{u, v) becomes a free fermion determinant 
(the sign counts the number of fermionic loops). There exists a well known 

formula in the theory of two dimensional free fermion models, relating T and J: 

(13) T(m, v) = ^ {J{-u, v) + J(m, -v) + J(-M, -v) - J{u, v)) . 

In the classical limit, equation J{u, v) = defines the spectral curve. Dynamical variables 
may be expressed in the terms of ^-functions on the Jacobian of the spectral curve. The 
sequence of Baecklund transformations, which is the "discrete time" in the classical model, 
is a sequence of linear shifts of a point on the Jacobian. The classical model was formulated 
and solved by I. G. Korepanov [3]. 

Classical integrability is based on an auxiliary linear problem. Equation J{u, f ) = is 
the condition of the existence of a solution of the linear problem. Our point is that in quan- 
tum g 7^ 1 case, the linear problem is still the basic concept of the solvability. Quantum 
J(u, v) is a well defined determinant of an operator-valued matrix, and J(u, v)|\l/) = is 
again the condition of the existence of a solution of a quantum linear problem. The poly- 
nomial structure of e.g. (f follows from a more detailed consideration of the quantum 
linear problem in a special basis of diagonal "quantum Baker-Akhiezer function" (related 
to a quantum separation of variables). 

The structure of the paper is the following. In the first section we recall briefly some 
basic notions of the classical model [3]: the linear problem, spectral curve and details of 
the combinatorial representation of the spectral curve. In the second section we repeat the 
definition of the quantum model and its integrability [1,2]. In particular, our definition 
of the spectral parameters differs from that of [1]. Quantum linear problem, derivation 
of (fTT|) and properties of various forms of (fT^ are given in the third section. The fourth 
section includes an example. 
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1. The Korepanov model 



We start with a short review of the integrable model of classical mechanics in discrete 
2 + 1 dimensional space-time [3] . The main purpose of this section is to recall the relation 
between Korepanov's linear problem, spectral determinant and partition function for free 
fermion model. Another aim is to fix several useful definition and notations. 

1.1. Linear problem. Consider a two dimensional lattice formed by the intersection 
of straight lines enumerated by the Greek letters. Let the vertices of the lattice are 
enumerated in some way. 

Consider a particular vertex with a number v formed by the intersection of lines a 
and /3, as it is shown in Fig. ^ It was mentioned in the Introduction, an auxiliary 
lattice is a section of three dimensional lattice, the vertices on the auxiliary lattice are 
equivalent to the edges of the three-dimensional one. In Fig. Q the dashed lines are the 
lines of auxiliary lattice, while the solid sprout from the vertex v is the edge of the three 
dimensional lattice. 

Let four free C-valued variables 

(14) = {a^,K,c^,d^) 

are associated with vertex v. In addition, let C-valued variables ipa and tpj3 are associated 
with the ingoing edges, and C-valued variables ip'a i^'/s associated with the outgoing 
edges, as it is shown in Fig. ^ (a certain orientation of auxiliary lines is implied). The 




Figure 1 . Vertex v is formed by intersection of a- and /3-lines of auxiliary 
lattice. Vertex linear problem is the pair of relations binding four edge 



variables. 



local linear problem is a pair of linear relations binding the edge variables. Its standard 
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form, the right hand side of Fig. ^ in matrix notations, is 



(15) f = X[AA { *■« ] , where A'l^.] S f V 

1.2. Korepanov's equation. Equations of motion in an integrable model arise as an 
associativity condition of its hnear problem. To derive their local form, consider a vertex 
of three dimensional lattice (not necessarily the cubic one), and sections of it by two 
auxiliary planes, as it is shown in Fig. |21 From three dimensional point of view, the 
auxiliary linear variables belong to the faces of 3d lattice, while the dynamical variables 
belong to the edges of the 3d lattice. Therefore are distinguished from A'^, but 
linear variables on outer edges ipa^ • • • ; "^75 iii both top and bottom auxiliary planes are 
identified. Consider the bottom plane first. The linear problem rule (|TH|l may be applied 




Figure 2. Left and right hand sides of Korepanov equation. 



three times for excluding internal edges, as the result one obtains an expression of the 
"primed" linear variables in the terms of "unprimed" : 



(16) 



where (cf. the ordering of 7 in column vectors) 
(17) 

( a2 62 \ 

-^0,7 [^2] = 



I ai hi \ 

Ci di 

V 1/ 



1 

\C2 d2 I 



/ 1 \ 

03 63 
V C3 j 



The top plane of Fig. El may be considered in the same way, 

where the matrices X^^^ are given by (fTTj) with A'^ = (a^, b'^, c^, d'^). 

The associativity condition of linear problems (jl6j) and (jl8|) is the Korepanov equation 



(19) X^AA] ■ XaAA] ■ X^rfi^s] = • X^,^[A',] ■ X^AA] , 

relating the set of 12 variables A^ with the set of 12 variables A'^, v = 1,2, 3. Equation 
(|19|) describes a single 3d vertex. Equations of motion for three dimensional integrable 
model is the collection of equations (jl9|) for all vertices of the 3d lattice. 

Korepanov equation needs a very important comment. Matrices X^^^[Ay] by definition 
(fT7j) act in the direct sum of one-dimensional vector spaces labelled by the indices a, /3, 7, 
etc. The matrix XQ,^/j[y4i] in the block (a,/3) coincides with ()15p. and in the block 

(7, . . . ) it is the unity matrix. In what follows, such "direct sum" imbedding of 2 x 2 
matrices X into higher dimensional unity matrices will always be implied. 

1.3. Linear problem with periodical boundary conditions. Korepanov's solution 
of the equations of motion is based on the solution of the linear problem for the whole 
auxiliary lattice. Consider the square lattice with the sizes x M. Let the lines of the 
lattice are enumerated by 

(20) an and pm , n = l,2,...,N, m = l,2, ...,M. 

A fragment of the auxiliary lattice is shown in Fig. |21 Notations for vertex and auxiliary 
variables for (n, m) vertex of the plane are shown in Fig. EJ The local auxiliary linear 
problem ()15|) for (n, m) vertex takes the form 

, (m-l) \ ( i{ra) 

(21) I \ =X\AnA ■ I ^;") 1 , n = l,...iV, m=l,...,M. 
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Figure 3. A fragment of the auxiliary square lattice. 
The linearity of the whole set of ()21|) with respect to ■i/'-s makes it possibly to define the 



^ ^t:'' 



An 1 



Figure 4. Notations for (n, m) vertex of the auxiliary lattice. 



quasi-periodical boundary conditions for them: 



(22) 



V^ir = and <„^^^=^;^, 



(n) 



where u and v are C-valued spectral parameters. 

Linear equations (pT|l may be iterated for the whole lattice as follows. Let 



(23) 



(m) 
a 



4?^ 



Then the repeated use of (PT|1 gives 



I I ) ^ ( ? 

where, in the terms of matrix imbedding discussed right after the (N+M) x (N + M) 
monodromy matrix lia.,i3 inay be written as 

(25) ^a,f3 = n n ^an,/3m[^n,m] , 

n m 

where 

(26) 11'^" ^ /l /2 • • • /aT-I /at , JJ^ /m = /l /2 • • • /m-1 /a/ • 

n m 

The boundary conditions (j^^ give = u^/'a'', i})^^^ = vil)^^\ so that becomes 

(^--- (::))(:?) 

The whole hnear problem has a solution if and only if 



(28) J(m, v) = det I 1 - X^..^ 



u 
V 



is zero. Equation J(m, t>) = defines the spectral curve for the model, equations of 
motion (fT^ for the whole three dimensional lattice have an exact solution in the terms 
of ^-functions on the Jacobian of the spectral curve [3]. 



1.4. Free fermion model. The determinant (jzHj) has the very well known combinatorial 
representation. Usual way to derive it is to define the determinant in the terms of the 
Grassmanian integration and then to turn from normal symbols to matrix elements. 

Let in this subsection ip and ip be the Grassmanian variables with the integration rules 
j dip = J dtp = and J ipdip = J ipdip = 1. Then the determinant PH|) may be written 



as 



(29) J{u, v) = I e^'^'^J Vijj'DiP 

where the "action" is 

N,M ( / j(m) 



n,m=l k \ r/3„ 
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and the measure is 
(31) 



n,m=l 

Spectral parameters appear in ()29|) via 



(32) 



-(0) -(M) -(0) -{N) 

i^a„ = and = vijp^ 



In the terms of Grassmanian variables, the exponent of a quadratic form is a normal 
symbol of some operator L, 



(33) 



exp|(^,,V';3) -XiA,] 

The fermionic coherent states are defined by 

(34) = |0) + |1)V', (^1 = (0|+V^(1| , 
and the extra summands in ()3U|) correspond to the unity operators 

(35) 1 = / 1^) e^^ # dip . 



It is important to note that the indices of the operator La,/3[Av] label copies of two- 
dimensional vector spaces (jMj) 3 x|0) + y\V) ■ Thus, L^^p acts in the tensor product 
of two-dimensional vector spaces, while Xa^p acts in the tensor sum of one-dimensional 
vector spaces. In the basis of the fermionic states 



(36) 



\na,rip) 



0,0),|1,0),|0,1),|1,1) 



operator La^p ()33j) is 4 x 4 matrix 



(37) 



-^a,/3[^v] 



1 














civ 


K 








Cv 


(iv 

















def 



where = Kc^ — a^d^ 



Besides, the Korepanov equation (fT^ is the equality of the exponents of the normal 
symbol form of the local Yang-Baxter equation 



(38) 



smce 



-^a,/3[^l]-^Q,7[^2]-^/3,7[^3] " -^/3,7 t^s] -^",7 [^2] -^",/3 [^3 



3J ' 



(39) 



\La,/3La^^Li3^^\ 



exp {V' ■ Xa^pX^^^Xp.^ ■ ijj} etc. 
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Turn now to the expression of the determinant ()28|1 in the terms of operators L. Let 
2N+M ^ 2^+^^ matrix L^,^^ be the ordered product of local L-s: 

(40) L«,/3 = nn La„,pJA^,m] . 

n m 

This is related to the monodromy matrix by means of (cf. (jSSl)) 

(41) exp I {ip^, ipf^) ■ Xc,,/3 ■ j I " ip(3) ■ 
Define now the boundary matrices for L^,^, 

(42) D{u) ^ ° ) , D^u) = n^".(^)' Df,{v) = llD^M, 
and let 

(43) T{u,v) = TTace(^Dc,{u)Di3{v)'Lc,,f3^ . 

By the construction, T{u,v) is the partition function for a free-fermion lattice model 
with the inhomogeneous Boltzmann weights - matrix elements of La„,f3^[An^rn] ^ and 
u, w-boundary conditions. It is the polynomial of u and v: 

N M 

(44) T{u,v) = 

n=0 m=0 

Sometimes a pure combinatorial representation of the partition function is very useful. 
Any monomial in T{u,v) PHjl corresponds to a non-self-intersecting path on the toroidal 
lattice. A path may go through a vertex in one of five different ways as it is shown in Fig. 
El (or do not go through at all). A factor is associated with each variant, these factors 
are the matrix elements of Lq^^[Av] (j?fj) . A monomial in T{u,v), corresponding to path 
C, is 

(45) tc= n f^- 

along path C 

Any non-self-intersecting path on the toroidal lattice has a homotopy class 

(46) w{C) = nA + mB, 

where A is the toroidal cycle along the a-lines, and B is the toroidal cycle along the 
/3-lines of Fig. El Then the element tn,m of (jl^ is 

(47) tn,m = J2 

C : w{C)=nA+mB 
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Ztj d-^j by 



Figure 5. Six variants of bypassing the vertex. Vertex factors /v are 
matrix elements of L. Note, in the variant the path is not self-intersecting. 

The determinant J{u, v) ()28|) is related to tn,m via 

N M 

(48) J{u,v) = 

n=0 m=0 

where the sign (^_)""^+"+'" counts the number of fermionic loops on the toroidal square 
lattice. The determinant J may be expressed in the terms of T and vice versa: 

J{u,v) = ^{T{-u,-v) +T{-u,v) +T{u,-v) -T{u,v)) , 

(49) 

T{u,v) = -{J{—u,—v) + J{—u,v) + J{u,—v) — J{u,v)) . 

The last equality is very well known in the two dimensional free-fermion model as the 
formula relating the lattice partition function and fermionic determinant. 

Now we may finish the collection of notions and definitions of the Korepanov model of 
classical integrable dynamics on three dimensional lattice and proceed to the description 
of their quantum analogues. 

2. Quantum model 

In the previous section, we did not pay any attention to the structure of vertex variables 
(fT^ - they were defined simply as the list of elements of X ()15|) . The key point for 
the quantization of the model is that it is possible do define a local Poisson structure on 
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[1] such that the transformation 
(50) A-^^A^^A^ A[ 0^2® A3 , 

defined by Korepanov equation (fTU|) . is a symplectic map. Symplectic structure admits 
an immediate quantization. We will skip here all the details and proceed directly to the 
Ansatz for quantized A^, X[A^] and //[A^]. The aim of this section is just to give precise 
definition ofT{u,v) ©. 

2.1. Quantum Korepanov and Tetrahedron equations. The local g-oscillator alge- 
bra 7i is defined by (0). The Fock space JF representation for g-oscillators corresponds to 



(51) 



Spec(h,) =0,1,2,... V 



Quantized dynamical variables (fT^ are the g-oscillator generators and a pair of C- valued 
parameters, ~ {'H^', Av,/iv): 



(52) 



Quantized X flT^ and L (jHTj) are given by 
(53) X[A] = 

and 



(54) 



^ 1 ^ 

Aw"^' yv 

-g-^Av/iv y 

One may verify directly, the quantum Korepanov equation 



v 



(55) 



is equivalent to the auxiliary tetrahedron equation (the quantum local Yang-Baxter equa- 
tion) 



(56) 
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Here the intertwining operator R = R123 acts in the tensor product of representation 
spaces J-'i®J-'2®J^z of three g-oscillators 7ii^2,3- Parameters Av, yUv of are the parameters 
of R. Both equations (|^ and (|^ are equivalent to the following set of six equations: 
(57) 

g^^xaya^ R , R x^ = f x^Xg + -^q^^+^^^^ R , 

^1/^3 



R g'^^xi 


A2 j 
" A^' 




Q 

Al/i3 


Rgh^xg 


- ^\ 




AiAi3 









g'^^YiXa j R , R = (y,y, + Ai/igg^^+'^^ya j R 



Classical equations (fT^ and (jHHj) follow from and ()56|) in the g — 1 limit of the 
well defined automorphism A'^ = R^^R^^. For irreducible representations of H^, R is 
defined uniquely, its matrix elements for the Fock space representation are given in [1]. 
Remarkably, Fig. |21 may be used for the graphical representation both of classical and 
quantum equations, in the quantum case the solid 3d cross in Fig. |21 stands for R. 

Integrable model of quantum mechanics may be formulated purely in the terms of 
matrices L fj54|) . The intertwiner R is related to evolution operators, it is another subject 
and we will not consider it here. Below we recall the definition of integrable model of 
quantum mechanics from [1,2]. 

2.2. Transfer matrix T. For the square lattice x M of the previous section, define 
the "monodromy" of quantized L-s literally by PUI): 

(58) L„,^ = nn Lo.„,pJ-^n,m] , 

n m 

and its trace (cf. (O) 

(59) T(M,t;) = 1ia.ce(^c,{u)Df3{v)l.c,^ , 



where boundary matrices D are defined by (jl2|) . To distinguish the classical and quantum 
cases, we use the boldface letters for the quantized T and its decomposition (jlH): 

N M 

(60) T{u,v) = 

ra=0 m=0 

Since the arguments of L-s are local g-oscillator generators, T(u, v) G 7i®^^ is by def- 
inition a layer-to-layer transfer matrix, its graphical representation is again the classical 
Fig. El The model is integrable since 

(61) T{u,v)T{u\v') = T{u,v')T{u,v), 
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(62) 



La,l3[-A.o] 



i.e. the coefficients t„ ^ in (jfiOj) form the set of the integrals of motion. 

Now we give the technical proof of the commutativity The commutativity of 

layer-to-layer transfer matrices follow from a proper tetrahedron equation [4]. In addition 
to La,f3[A^] dSH), define 

^ 1 \ 

Ao(-g)^° yo 
g-^o/UoXo fJ^o{-q)^° 
y g~^Ao/io J 

where ~ ('^o! Aq, /Uq). The constant tetrahedron equation for L and L, 

(63) La,a'[-Ao]L/3,/3'[Ao]La,f3[A]La',/3'[A] = Lq,',/?' ['^]-^/3,/3' [^o]-^Q,a' [-^o] 

may be verified directly, it is just 16 x 16 matrix equation with the operator- valued entries. 
Another technical relation is 



La,a' [Aq 



ho 



(64) DMDa'{u)Lo.,^>[Ao] = 
where D is given by (021) • Combining for the whole lattice, we come to 

(65) '^a',13' = '^a',13' ^cx,l3 ^13,13' '^a,a' , 

where besides the "monodromies" (PUI of L[y4v] 

n,m n,m 

we used 

n m 

Multiplying (jU^ by Da{u)D(i{v)Dct'{u')D^i{v') and taking into account (jMjl . we get 



/ / \ ho 



) L^,/3' ■ £'«(m)^/3(^^)L«,^ ■ /^«'(m')^/3'(^')W,/3' 



(6J 



Dc.'(^i')-D/3'(^')W,/3' ■ -Da(M)-D/3(^^)ILa,/3 



ho 



t> \ ho 



Now taking the trace over the representation space JFq of Tio and denoting 



(69) 



^a,l3,OL' ,13' 



Trace ( Lq, 

J^o V ' 



■u V 
uv' 



ho 



L 



16 S. SERGEEV 

we come to the final similarity relation 
(70) 

and therefore two transfer matrices T(u,v) = Trace! Da{u)Di3{v)'La (3 ] commute. 
In the limit g — 1 coefficients of (jUUI) become the involutive moduli of the spectral curve 



). Since the moduli of the classical spectral curve are independent and the spectral 
curve defines completely the solution of the classical model [3], the set of integrals of 
motion t„ ^ is complete. 

The combinatorial representation for T{u,v) (j59|l is equivalent to the combinatorial 
representation for T{u, v) (jlHjl . The vertex factors of L[A^] (jHTjl in Fig. Elare to be replaced 
by corresponding elements of L[A^] (fS^ . Equations and (jTTj) remain unchanged. In 
particular, with the help of the combinatorial representation one may easily see 

(71) T{u,o) = n (i+^n^".-^''""") ' '^(0'^) = n (i+^ri/'".-^''"") ' 

n \ m / m \ n / 

SO that the basic integrals of motion are q'^" and g'^™, 

(72) J7n = ^ ] ^n,m ^ud /C^ = ^ ] ^n,m ■ 

m n 

Eigenvalues of Jn and Km fix a subspace in the state space of the model. 



2.3. Transfer matrix (j59p and 2d quantum inverse scattering method. The trans- 
fer matrix may be identically rewritten as the trace of 2d monodromy matrix 



(73) T(M,t;) = Tia.ce [Df^iv) 

^ V n 

where 

(74) 2.f{u) = TT^ce(n^^{u)\[K^,pJA.,rn]\ 
Equations ()fi3j) and ()fi4|l provide in particular 



ho 



f(o) . \ 7(0) . n r,,u(™) . n r?/u(™) 



(75) 



/\ ho 

U 
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and therefore 
where 

(77) ^^^A^/u') = Trace H^") "HV^/^jJA] 

L ^ ^ m 

In the same way, the other direction of the lattice may be chosen, and the dual Lax 
operator 

(78) ^i^\v) = TTs.ce(D^Jv)flL^„,pJAn,rn]] 
may be considered. 

Further in this subsection we recall the structure of 2d Lax operator (f7^ and its R- 
matrix (f77j) . Matrix , as the matrix in the auxiliary space = (c2)®A'-f ^ith operator- 
valued matrix elements from 7i®^^, has a block-diagonal structure. Let, in the concordance 
with and (jHBj) . |0) and |1) be the basis of C^. Define ipo, ipj, fjj', etc. as the following 
elements of V^: 

= |0)®---®|0) , 

(79) 

(Pj = |0)®--- |1) ®---(8)|0) , 

j place 

etc., in general 'Pj-^,...j,^ has |1) on jf", ... , places, ji < ... < jm- In the basis of (p 
(80) 4")(m) Po = V>o , where Lo{u) = 1 + uY[Xn,^q''"''- . 

m 

Next 



11) si^;\u)pk = 



where Lj^k{u) are matrix elements of Lax operator for the vector representation of ^^(sZm), 
they are given in [1]. In general. 



M 

[82) £b(u) = e L^^(u) 



m=0 



where block Um corresponds to jn\{M-my. -dimensional vector space tt^^^ G V/3 with the 
basis <^ji,...,j„, ji < ■ ■ ■ < jm- 
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Matrix 9^^,/3' (fTTjl has the block structure as well, 



(83) 



M 



m,m'=0 



where Ao,/io are extra parameters of dZZj), and Ruj^^ui^, is the Uq{slM) -R-matrix for 
the representations n^)^ ® tt^^^, • In particular, in the sector vr^-^ ® vr^^ with the basis ipj 
(jZni), one can obtain the fundamental i?-matrix (we used the Fock space (jHT|l for Tio for 
the calculation of the trace in (fTTjl ): 



(84) 

where 
(85) 



{u-q'){u-l) ^ ^ ^''^ 



1, /?y = g-l(^i-g2 



Rj^ j{u) = q ^(1 — g^) for j < k and i?; 



q ^u{l — (f) for j > k 



(m) , 



Lax operator (u) (j74p has the center g^" (j72j) . while g'^'" is the center of £,', 
(fTHj). The quantum space of £^ is J^®^^ it may be decomposed as 

oo 

(86) = e TTj^, , 

J=0 

where J is the eigenvalue of JT, and ttj^^^ is the rank- J symmetrical tensor representation 
olUq{slM) (dominant weight Jui). 

We would like to conclude this subsection by the example of M = 2 containing the 
six-vertex model. The block-diagonal structure of £ ()74|) is 

/ 1 + uAiA2g''^+^2 ^ 

(87) £(«) = L{u) 

V /ii/i2(g'^^+'^^+g-'MAiA2) / 

where index n of (f?^ is omitted, and 2x2 central block is 

^.i{q^^ - u\iX2q^^'^) -g"^MAi/iiXiy2 
-g"ViA2/U2X2 /i2(g''2 - MAiA2g'^^"^' 

Fixed integer J7 = hi + h2 in the quantum space corresponds to spin J /2 representation 
of 5/2- For spin 1/2 representation {J = 1), matrix elements of (|HH)) may be presented by 



(88) 



L{u) 



(89) 




h2 




xiy2 





l-g2 



, yiX2 
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and in the homogeneous case Ai = A2 = /Ui = /i2 = 1 the matrix ()88|) becomes exactly the 
six-vertex i?-matrix. 

Let the chain of Lax matrices ()74|) is given, n = 1, 2, A^. Six vertex model corresponds 
to the choice Jn = ^ for all ^- The values of two extra integrals of motion, /Ci and /C2 
of (f?^ . are related to total spin of the chain: /Ci is the number of "spins up", /C2 is the 
number of "spins down". 

3. Quantum curve 

Solution of the classical equations of motion is based on the notion of the spectral curve 
J{u,v) = (I2HI)- In the previous section we succeed in construction of the "quantum 
partition function" T{u, v) producing the set of the integrals of motion, but we did not 
answer the question: what is the quantum analogue of J{u, v) = 0? 

The answer is the following (we are repeating the Introduction). Let u, v be an addi- 
tional auxiliary Weyl pair, 

(90) W : u V = g2 V u 
serving two variants of "quantum-mechanical" notations 

{Q\u) = Q{u) , u\u) = \u)u , vIm) = \q^u)v , 

(91) or 

{v\Q) = Q{v) , {v\v = v{v\ , (t>|u = u{q'^v\ . 

For the given layer-to-layer transfer matrix T{u,v) ()59I60|) . define 

N M 

(92) J(u, v) = Y.ll i-qy-'i-un-^rr tn,m . 

n=0 m=0 

One may easily see, in the limit q ^ 1 (jU^ becomes ()48p. Operator J(u, v) belongs to 
j^mM ^ yy;_ | ^ jr^M eigenvector of T(u, v), 

(93) t„ „ I I tn,m ; 

SO that 

(94) J(u,v) ( |t) ®ic/ ) = ( |t) ® irf ) J(u,v) , 



where J(u, v) G W, fjll|) . Define a linear space \1' by 

(95) J(u,v)|^) = or (^|J(u,v) = 0. 
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The linear space may be decomposed with respect to the basis 1 1), 



(96) 



1^) 



t)®|Q,) or (^1 



{t\ ® {Qt 



The point is that {Qt\ and \Q^) have a simple and predictable structure, in particular 
{Qt\u) and {v\Q^) may be defined as polynomials. 

Functions {Q\u) and {v\Q) (index t is usually omitted) obey the linear difference equa- 
tions dnn), 



(97) 



(^1 J(u,v) I t) O |m) = and (t| O J(u, v) |^) = 



which are exactly Eqs. ((Tj) and (jH)). 

In this section, we will derive (I92j) and (j95|) as a solution of quantized linear problem 
P5|) for the whole auxiliary lattice. A non-trivial problem is how to define a quantized 
linear problem in such a way, that each linear variable belongs to the same subspace 
(subspace (j^^j) in the final solution). This is the problem of consistency of linear problem 
and invariance of quantum curve. 



3.1. Generalized matrix of coefRcients of linear problem. Consider the linear prob- 
lem fllSp on an arbitrary lattice with u, f-periodical boundary conditions. The set of local 
linear problems 



(98) I ' 1 = X. ■ I ' " 1 , X. 

may be rewritten for the whole lattice in a matrix form, 

(99) E^^.^-^^- = 0- 



(100) 



3), corresponding 


to dHHl 


, are 


/ .. : 










... \ 


1 





— Ctv 


-K 








1 


-Cv 







... 





1 







... 








1 














■•• / 



for %l) 



( : \ 

V ''■ ) 
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In the quantum world, matrix elements ^k,j are non-commutative operators ikj- Therefore 
fl99p may have two slightly different forms, 

(101) A). 5^4,il^i)=0 or B). ^(^j|4,i = 0. 
Let the operators ik^j obey the following exchange relations: 

^k,j^k' ,j' ^k' ,j£k,j' (-k' ,j'(-k,j ^kj'^-k' ,j ; k ^ k j j 7^ j i 

(102) and 

(■k,j(-k',j ^k',j(-k,j • 

The aim of this subsection is to establish fundamental properties of such i and the form 
of solution of (ITTITI A) and B). 
Define 



(103) det£ J2(-rl[ia„, = $^(-)X„i^..,2C3,3 ••• 

where a are the permutations of the indices 1, 2, 3, . . . . According to the first relation of 
fll02|) . the definition ()103p is invariant with respect to the ordering of j, for instance 

(104) ^(—)'^d, 1^2,2^3,3 ■ ■ ■ = ^(—)'^C2, 2d, 2^3,3 ■■■ = •••, 
and in general 

(105) det€ = E(-)'^n^-.'-. 

o- j 

where r is any permutation of 1, 2, 3, ... . Define next the algebraic supplements Aj ^ of 
£ (adjoint matrix) as ()103|l -determinants of the minors of £, 

(106) ^Aj^k£k,j = ^ik,j Aj^k = dete. 

k k 

The first equality here follows from ()105|) . the last one in the definition of Aj^. The 
second line of fll02|l provides 

(107) ^Aj-fc4,i' = 5^4j'A,-fc = if jV/- 

k k 

Therefore (det €)^"'^Aj ^ and Aj ^(det £)^-'^ are two variants of inverse matrices, 

(108) J2{dete)-'Aj,k£kj' = 5^4yA,-fc(det£)-^ = Sjj, , 
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or, since the inverse matrix must be both left- and right-inverse, 

(109) 5^4j(det€)-iA,-fc„ = ^A,-fe„(det£)-i4,i = 4,A:o • 

j j 

The main property of the elements of inverse matrices is the commutativity of their matrix 
elements with the same ko, e.g. for the variant A), 

(110) [(det^)-iA,-fc„ , (det^)-iA,,fc„] = 0. 
It is a particular case of the following statement: 

(111) ^ 4j rrij = Ek , SkeC =^ [mj,mj>] = 0. 

j 

To prove (jlllll . consider 

(112) Ck,k' = ^(4,j4',j' - £k',j£k,j')(mjmy - my-nij) , k < k' . 

j<j' 

Due to the first relation of ()102|) and definition ()111|1 . Ck,k' = £k£k' — £k'£k = 0. From the 

(2) 

other side, let A^. ^.,|^^, is the matrix of the second algebraical supplements of £: 

(113) ^f)\kA^KA' j' - ^k',A,j') = Sij6i,j> det£, i<i' , j < f . 

k<k' 

Then = ^flik k''^k,k' = det£ ■ {mjUij' — mj/rrij), what proves the commutativity 

k<k' 

of nij. The commutativity (jllO|) corresponds to Ek = 5k,ko- 
For Ek = 5fc,fco with fixed /cq; let 

(114) m^.^.^ = (m.-J-^m, = A]lk,-^,,k, , = ^.mAJoM ■ 

Then the set mjj^ is commutative (/cq is fixed), as well the set mj-jg, and besides 

(115) det£ ■ m^.^.^ = m^-^-^ • det£ . 
Inverse relations p09p give 

(116) ^4,im^j^ = 5k,koA-]ko det-^, ^m^.^.^4j = 5k,ko det i A']^^ . 

j j 

Equations (|116|) give the formal solution of the linear problem (jl01|) A) and B): let {il^jo) 

or {ipjgl for some jo are defined by det£ ■ li'jo) = or {ipj^l ■ det£ = 0. Then 

(117) A). liPj) = mjjjV^j,) or B). {iPj\ = {ipj,\m'jj, ■ 
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Common feature of all \ipj) or {ipj\ is 

(118) A). det£ ■ = or B). ■ det£ = V j. 

3.2. Extended algebra of observables. To establish the relation between the algebra 
of matrix elements ()1()2|1 and our case of X ()53j) . we have to introduce the extended 
algebra of observables 2tv = Tiv ® VVv. Define 



:ii9) 



hv 



where g'^'', x^, yv are generators of g-oscillator Ti^ ((T)) and X^, fi^, are the generators of 
Weyl algebra W^: 



(120) 



The g-oscillators and the Weyl algebra elements for different vertices always commute. 

Matrix £, defined according to ()10U|) for vertex matrix X[2lv] pi9|) . belongs exactly to 
the class ()102|) since 



:i2ii 



i^x uy — nq^ = wy i/x — Xq^ , 

and 

[uy , Mg"^] = [i^x , Ag*^] = 

and elements of 2lv for different vertices commute. Determinant of £ has a combinatorial 
representation of Fig. eqs. ()45f47p48j) with 

/ 1 \ 
Aw*^^ i^vYv 



(122) 



y 

Let 

(123) det€ = J[2t] = ^(-)"™+"+'"M"t;'"J„,„[2l] , 

where u and f are C- valued spectral parameters introduced according to (j22|l . Let locally 

(124) 

Av = A^e^^ = fi^e^" and therefore ul = -q'^X^fi^e^''+^'' , [Q^, P^] = logq^ . 
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In the combinatorial representation of the determinant (Fig. eqs. ()45p47f48|l ). consider 
a path Cn,m of the homotopy class nA + mB. A monomial summand, corresponding to 
this path, may be factorized as 

(125) Jc = tc e'^^^"'") . 

Monomial t^^^ gathers all the g-oscillators and C-valued parameters Av, fi^ and — g^^Av/iv, 
and therefore it is exactly the C„ ^-monomial of T{u, v) (up to unessential renormalization 
of Xv and yv). Operator is a sum of local and P^. One may easily see, 

(126) [0(C„,J,0(C;,^,)] = inm'-mn')\ogq' y C,C' . 
Let Cifi and Co,i be two particular fixed paths. Due to p2(j|l . all 

(127) 4>{Cn,m) = (f>{Cn,m) ' n4>{Cr^,) - m0(Co,i) 

commute. Therefore one may diagonalize them simultaneously and without lost of gen- 
erality (since Av, /iv are free) put 

(128) 0(a,™) = . 
On this condition, 

(129) 4>{Cn,m) n0(Ci,o) +m0(Co,i) = nQo + mPo, 

and exponent of 4>{Cn,m), together with the combinatorial factor (_)"'"+"+'"^"^™-^ be- 
comes 



(130) (^^Y'^+n+m QuQo+mPo ^n^,m ^ -nm ^ ^ _ gQo^ , V = 6 



Po 



V . 



Therefore, on the subspace (|128|) . J[2t] becomes exactly J(u, v) (jU^ . 

Definition of u, v p30|) on the subspace p28|) may be written in the terms of Av, fi^ as 



(131) M]^A„r« u]^A„r«, ^nMnm'^" V J]^ 



m TO 
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3.3. Structure of a solution of linear problem. Turn now to the structure of linear 
spaces and defined by 

(132) A). J[2l] 1^) = and B). (^| J[2l] = . 

According to pi()|114|115|l . each of and belong to these subspaces, and there 
exist sets of commutative operators rrij and nij (I117p such that \ipj) = nij |\1') and 
{ipjl = (\[^| m^-. In particular, one may consider the linear spaces and (\E'| in the basis 
of diagonal rrij , : 

(133) A). = |^)mj , or B). {ipj\ = mj{<i!\ , 

where ruj are eigenvalues of rrij etc. The every local pair of linear equations (j98|l take the 
form 

(134) ^ 

A). { or B). I 

[1 - i/vx,a - M,g^^Dv)|^) = [ (^1(1 - a^^vXv - D.fx.q''^) = 

where in the notations of (jnH|l = — —, = , = — — and = . 

rria' nia' rnpt mpi 

The hnear spaces |^) and (^| belong to a module of U®^^ ® y\;®NM _ "physical" 
part is jF®^^^, the artificial Weyl part may be defined in many inequivalent ways. Let \n) 
be the state of jF®^^^ with bosons in the vertex v, n = {n^}. We will focus on the 
formal states of W^^*^ for (IT!?^ in the basis (HHSI): 



(135) A). = y{n\®idj\^!) or = {m\\\n) ^ id] . 

Equations ()134p give for p35|) 

(1 - AvAv)|^o) = Vv and 



A). <; -rT ^v-(Mv^v;^^)n 



(136) or 



(^o| (1 - Atv^v) = V V and 



Here we used the representation 



(137) X |n) = |ra-l)Vl-g2n^ y\n) = \n^l)^l-q 



2+2n 



26 S. SERGEEV 

Turn at the first to the conditions for \l/o in (jlH^j) . Since the total Fock vacuum is the 
eigenstate of g-oscillator counterpart of J [21], one may consider directly 



A) . = \^{0\^idj J[2l] 1^) = JoiW®"^*'] l^o) 
(138) or 

B) . = (^|J[2t] (\0)(S)id] = {^o\ 3o[^^^^'] 



where 

(139) Jo[w«n = 11(1 - ^n^-) 11(1 -^n^^ 



nmj 

m 



Evidently, Jq commutes with all (f) p27p and therefore allows the projection W®^*^ W 
(IT2H1)- Simply applying 

(140) Jo[W®^"Vh^° Jo(u,v) = \{{l-^\{^^nm)\{{l-^^\{Km) 

m n n m 

Now the state of the Weyl algebra u, v may be chosen in the most convenient way. For 
the case A the proper basis is {v\ of (PT|) . and for the case B the proper basis is \u). If 
u = Un for the case A or t> = f ^ for the case B, 

(141) Un = Y[\~^, Vm = Y[fi~^, 

m n 

then 

(142) A). (t;|^o) = 1 or B). {^o\u) = 1 . 

Let further \t) and {t\ be the eigenstates of tnm, (|^- Analogously to ()135p let 

(143) A). I^t) = ({t\®id]\'^) or B). = (^|(|t)®id 



It follows from the second lines of 

(144) A). 1^,) = Vtiu,^)\^o) or B). (^,| = (^oin'('^>A), 

where e.g. Vt is a polynomial of fi^ of a total power not more that the number of bosons 
in its structure with respect to is rather simple. In addition, since we consider the 
eigenstates in the Fock counterpart of 2t®^^, polynomials Vt and must commute with 
all (fnTjl and therefore must allow the projection VV®^^'' W (fT^ : 

(145) Vt w-^Pt(v) and V't P;(u)w-^ , 
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where = — vu, w-factor comes from factors; J and K are integers; Pt and P/ are 
polynomials of a power not higher than the total number of bosons in the state 1 1). 

Taking now into account ()142|) . we come to the final statement: in the restricted algebra 
y^^NM ^ yy ^^iQ cquatious 



(146) A). ® (t;|J J(u,v) 1^) = and B). (^| J(u, v) O |u) J =0. 
have the solutions 

A) , li u = Un then = f"'^"/^Q„(f) 

(147) and 

B) . ii V = Vm then = u'^'^^'^Qmiu) 

where Un,Vm, are given by ()14H) . and the integers J„, Km and degrees of polynomials Qn, 
Qm are not higher that the total number of bosons. 

In the next section, considering the examples, we will see that J„ and are the 
eigenvalues of J'n and )Cm (|7^ and, moreover, the degrees of and Qn are exactly 
and Jn- 

In the next section we will add and g^'^™ to the definition of m„ and v^- This 
allows one to cancel the half- integer pre-factors ti^-^"/^ of Qn{v) and u'^"^^"^ of Qm{u). 
The corresponding values of and Vm may be obtained via conditions Q{0) = 1 and 

Qn(0) = 1. 

4. Examples 

Let us illustrate ()97j) for the six- vertex model first, and then for an arbitrary square 
lattice. 

4.1. Six- vertex chain. Consider the lattice with M = 2 and arbitrary A^. This is the 
case of six- vertex model, where (jU7|) becomes the Baxter equation. It follows from ()87|). 

n n 

(148) XI = ^^^^ ' 

n 
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where t{u) is the transfer matrix for the Lax operator ()88|) . Its opposite elements are 

to,i = n {^'n^q''"■') + n (/^n-sg*^-^) , 

n n 

(149) 

n n 

Applying the rule (— u)"'(— v)™, we come to 

(150) 

J(gu, v)v~^ = — t(u) 

+v-^ 11(1- q-'nXn,iXn,2q'"' ''''"' ') + v J] /x„,i/i„,2 (g^-^+^-^ - q-'uXn,iXn,2) . 

n n 

Equation (Q| J(gu, v)v^^|n) = is exactly Baxter's equation for Uq{sl2), 
(151) Qiu)tiu) = Q{q-\)(t){u) + Q{q\)(t)'{u) 

where 

0(m) = v-^\[{l-q-\Xn,iXn,2q''"^'^''"^') , 

n 

(152) 

and 



b'(u) = V 



Yl^^n,l^^n,2 (gh"^l+^-2 - q~^uXn,lXn,2) • 
n 

Condition Q(0) = 1 gives (see (ITm i 

(153) t'-' + t^n/^^'i/^-'^g''"''''''"'' = Ilf'n,iq''"'' + H'^^'^g'^"'' ' 

n n n 

which has two solutions corresponding to two Baxter's functions Q: the first one is 

(154) v = v^= (llfin,iq''" 



so that Q = Qi{u) is a polynomial of the degree h„^i = /Ci (see fl7^). this value of 

n 

the degree follows from the second line of eq. 1)149^ : the second solution is 



1 

def 



(155) V = v,"^ \\[^^^,,q''-^- 



n 



and corresponding Q = Q2{u) is the polynomial of the power ^^hn,2 = A^2- In some 

n 

sense, two functions Q correspond to two sheets of classical q = 1 spectral hyperelliptic 
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curve v~^(f){u) +v(j)'{u) = t{u). Note, we consider now the Bethe Ansatz equations for 
Uq{sl2) chain with arbitrary J7„, this is the inhomogeneity of highest spin. Six-vertex case 
corresponds to Jn = ^ for all n. 

In the inhomogeneous case the spectrum of t{u) follows from ()15H1 and just the condition 
Q ^ (without fixing v and considering the polynomial structure of Q). Let us fix 
J'n = I, the six- vertex case and a priori the Fock space representation. Simply substituting 
u = (A„,iA„,2)~^ and u = q'^{\n,i^n,2)~^ into (jl^lj) and excluding Q, one comes to 

with n = 1, . . . ,N. These equations and conditions ()149j) do produce the whole spec- 
trum of t{u) for inhomogeneous six- vertex model. 

Equations p56|) may be obtained in the other way - via Q{v) = {v\Q). Condition 
Q{0) = 1 fixes u = Un, 

(157) Ur, = (A„,lAn,2g'^-^+'^"'^)"' , 

corresponding Q = Qn{v) is a polynomial of the power J'n = h„ i + h„ 2- In the six- vertex 
case J'n = I, therefore Qn{v) = v — (n, and equating to zero each f-term of the right of 
equations (|97j) . one comes to 

JJ(1 - q^Un/Uk) 

(158) C = = - ' , , . 

k 

The second equality gives p56p . 

4.2. Uq{slM) equations. Let the lattice has arbitrary N and M. We will consider fj97p 
for {Q\u) = Q{u) and for {u\Q) = Q{u) - the last function was not mentioned before, but 
it is interesting to discuss what it is. 

Equations {Q\J{u,\)\u) = and (m| J(u, v)|Q) = read correspondingly 

M M 

(159) ^Q(g2™u)(-t;)"^r„(M) = and J] ^(^~'"")(-^')"'^(^~""^) = 0' 

m=0 m=0 

where 

N 

(160) Tm{n) = $^(-g)""^(-w)X™ , 

n=0 
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cf. ()5|fij) . Combinatorially, one can obtain the following summation formulae: analogue 
of 



n n ^ m ' 

(161) 

n n ^ m m ' 

and analogue of ()149j] 

m m ^ n 

(162) 

m, m n " ^ 



-1 



Let us fix notations for the inhomogeneities (cf. (^SJ) 
It follows from ()16ip . 

(164) ro(n) = n fi - -) ' ^*^('") = n^™ n - 

„ \ y m n ^ ' 



For the normalization (5(0) = (5(0) = 1 both equations ()159|1 are equivalent to the first 
line of (Una): 

M . . 

(165) J^(_^)-r„(0) = n 1-- =0 

m=0 m ^ ™^ 

Suppose, Q{u) ~ -u^ and (5(m) ~ when m ^ cxd. Then (I159p and the second line of 
P62|) provides the following conditions for K and K: 

V 



(166) JJ(l_g2M„_|^J 

and 

(167) n ( 1 - g-'^-'^-'^-— ) = . 



The charges J'n and /C^ are given by ((721). Equation ()165|) has the solutions v = Vm 
m = 1,2,3, ... , M. In what follows, we will assume the generic set of /i„m, so that all f„ 
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are different. Let v = Vm corresponds to Q{u) = Qmiu) and Q{u) = Q^{u). Then fll66p 
and ()167|) defines the leading u ^ oo asymptotic 

(168) Qm{u) ~ m'^™ and Q^{u) ~ m^™ , 
where 

(169) Km = -J2i^+Km) ■ 

n 

Thus, for the Fock space representation all Qm{u) are polynomials of the power JCm- 
Functions Q^{u) are rational functions, below we construct them in the terms of Qm{u). 

Instead of the Fock space representation Spec(h) = 0, 1, 2, oo, one may consider the 
anti-Fock space, Spec(h) = —1, —2, —3, — oo. Then fll34|) are to be solved in different 
way, and as the result Qj^{u) become polynomials and Qm{u) become rational functions. 

The dual "T-Q" equations for sl^ correspond to the evident exchange M, u v 
and J IC. 

Turn now to the form of Q^. for the Fock space representation. The detailed "arith- 
metical" consideration of ()159|) as a set of linear equations allows one to conclude for 
instance 

(170) QmW = ' 

where 



(171) Wuiu) = det\\Q^iq'''u)vi\ 



\i,j=l,...,M-l ) 



and ^ , = V1V2...VM — \ n ^ ■ AH the other O correspond to (11701) and (117111 with 
permuted set of indices of Qj. Inhomogeneity of Vm is important in this consideration 
since if Qm(0) = 1, then WMiO) = vi...vm-i Y\, ("^^ ""^i)' therefore WMiu) is not 

l<i<j<M 

zero. As well, for the generic Vm all Qm{u) are functionally independent, their Wronskian 



(172) W{u) = det\\Q^{q^^'-^\)vY%j=^_M 
obeys 

(173) = V1V2...VM- 



W{q^u) "' Tq{u) 

with the initial condition W{Q) = Y\ ("^j ~ "^'i); therefore W{u) 7^ 0. 

l<i<j<M 
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Using these "arithmetical" considerations, one may express the fundamental transfer 
matrices Tmiu) of sIm via determinants det | |<5j(5'^^'tt)f |ij=i,...,M5 where Pi is a subset 
of (0, 1, . . . .M). Moreover, for a more general sets of pj, any transfer matrix of sIm may 
be expressed as such a determinant [5]. 

The nested Bethe Ansatz equations (see [6] and e.g. [7]) may be derived from the 
generalized "T-Q" equations in an "arithmetical" way as well. Let 

^^^^^ Qraiu) = det\\Qj{q^'u)v]\\ij=i^...^m , 

Vm{u) = det\\Qj{q'^'u)Vj\\j=i^,„^rn;i=0,2,...,m , 

with Qo = I and Vq = 0. Function Qm{u) is a polynomial of the power /Ci + /C2 + - • ■ + /Cm- 
Equations 

(175) Vm{u)Qm^l{u) = V.m^l{u)QM + —Q,n~l{q^u)Qm{q~^u) 

for m = 1, 2, . . . , M — 1 are just the determinant identities for ()174|) . For m = M the 
system (I175jl should be completed by an equation, following from ()159|) : 

(176) Ti{u)Qm-i{u) = Tq{u)Vm-i{u) + vi- ■ ■VM-iTM{u)QM~i{q^u) , 

where T o and T m are given by p(j4j) . The nested Bethe Ansatz equations are a closed 
system of algebraic equations for roots of Qm{u) following from ()175|17(i|l . 

5. Conclusion 

This paper has a modest aim just a to give a correct form of "T-Q" equations. We can 
say nothing about their solution. But we would like to note, that from the point of view 
of three dimensional models, the thermodynamical limit is the limit N, M ^ oo with 
non-singular ratio : M. The nested Bethe Ansatz equations was never investigated in 
this limit since Qm{u) p74|) with finite m has a finite number of roots. In addition, an 
excitation corresponds to a change of the structure of occupation numbers. From sIm 
point of view of the previous section, it corresponds not only to a change of /C^ related to 
the powers of Q-operators, but as well it corresponds to a change of J'n which is a change 
of the s/Af-structure of the nested Bethe Ansatz. 

Let us better conclude this paper by a brief comparison of two exactly integrable models 
in 2 + 1 dimensional space-time. From the point of view of the algebra of observables, 
these models should be called "g-oscillator model" and "Weyl-algebra model". The last 
one is a quantum-mechanical reformulation of Zamolodchikov-Bazhanov-Baxter model of 
statistical mechanics, which has a long history [8-12]. Both models are based on two 
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slightly different forms of local linear problem, the linear problem for Weyl-algebra model 
may be found in any of Refs. [13]. Solution of both classical models may be expressed in 
terms of algebraic geometry, [3] and [14]. Equations of motion may be understood as a 
canonical mapping conserving certain symplectic structure, [1] and [15]. Poisson structure 
allows an immediate quantization, [1] and [13]. Quantum-mechanical integrals of motion 
may be combined into a direct sum of transfer matrices for fundamental representations 
of either sIn or sIm, [1,2] and [16]. And finally, the solvability of the models is based 
on quantized auxiliary linear problem and remarkable features of their operator-valued 
matrices of coefficients, this paper and [17]. 

The continuous limit of both classical models may be illustrative. Equations of motion 
for six fields q^, q*, j = 1, 2, 3, follow from the action 



where for the classical continuous limit of g-oscillator model V = qlq^q'l —■ ?iQ'2?35 this is 
nothing but the model of three- wave resonant interaction [18]. For the classical continuous 
limit of Weyl-algebra model the potential is = {ql — q2)iq2 ~ QsJils ~ Qi) [19]- 
Acknowledgements The author should like to thank V. Bazhanov, V. Mangazeev, M. 
Batchelor, M. Bortz and all the Mathematical Physics group of the Department of The- 
oretical Physics of RSPhysSE for fruitful discussions. 
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